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In this paper, load frequency regulator based on linear quadratic Gaussian 
(LQG) is designed for the multi-area power system (MAPS) with 
communication delays. The communication delay is considered to denote 
the small time delay in a local control area of a wide-area power system. The 
system is modeled in the state space with inclusion of the delay state matrix 


parameters. Since some state variables are difficult to measure in a real 


modern multi-area power system, Kalman filter is used to estimate the 
Keywords: unmeasured variables. In addition, the controller with the optimal feedback 
gain reduces the frequency spikes to zero and keeps the system stable. 
Lyapunov function based on the linear matrix inequality (LMI) technique is 
used to re-assure the asymptotically stability and the convergence of the 
estimator error. The designed LQG is simulated in a two area connected 
power network with considerable time delay. The result from the simulations 
indicates that the controller performed with expectation in terms of damping 
the frequency fluctuations and area control errors. It also solved the 
limitation of other controllers which need to measure all the system state 
variables. 
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1. INTRODUCTION 

Automatic generation control (AGC) [i.e., automatic voltage regulator (AVR) and load frequency 
control (LFC)] for the multi-area power system (MAPS) has made the system highly economical and reliable 
to generate and transmit active power to load at safety frequency level and nearly nominal voltage in the 
modern MAPS under deregulated environment. From studies, the dependability of the active power to the 
load depends on the power regulation frequency. Therefore, frequency regulator is very important to keep 
frequency at the safe range in the MAPS. The job of frequency regulator otherwise known as the LFC is to 
retain the regulation frequency over randomly active power load changes sometimes known as disturbances. 
Another duty of the LFC is to control the tie-line power exchange error [1] that might exist in the MAPS. For 
simplicity, multi-area generating units are connected through tie-line communication link in other to ease the 
LFC. The use of this link introduces a new error called area control error (ACE) into the LFC. For example, 
if active power loads change affects one area, because the wide area power systems are interconnected 
through communication link. Therefore, this will not only affect the traditional governor speed control of the 
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generating area but also create control problems in the other areas. However, it is necessary for the area 
subjected to load change control and balance the fault within its area without other areas being affected. 
Otherwise, there is going to be an economic and reliability conflicts between different power generating 
companies tie together for example, GENCOS 1 and GENCOS 2 generating areas interconnected. 
Hypothetically, separated load frequency regulators are needed in each area to act swiftly during faults and 
able to set their setting parameters. Moreover, ACE signal receive by each controller is through the 
communication link, respectively. Therefore, the goal of control engineers is to develop the LFC to damping 
frequency error and ACE to the safe point after the disturbance has acted upon on any areas. The old 
proportional integral (PI), proportional integral derivative (PID) tuning based the LFC scheme [1] were 
applied to regulate frequency in the MAPS in the early days. However, with the complexity in modern power 
system, these controllers were seen degraded [1]. Recently, the development was made to use some 
algorithms combined with PID for the complex power system which were discussed in [2]-[9]. As power 
system becomes more complex, tuning parameters for the LFC becoming more complication. Therefore, the 
controller design needs to be robust enough to MAPS disturbances. Guo [10] the author applied sliding mode 
control (SMC) to control frequency in complex power system. However, SMCs were having chattering effect 
as a result of some neglected tuning parameters in the system model. The uses of SMC combined with some 
techniques to decrease chattering and minimize frequency error/ACE were discussed in [11]-[20]. Linear 
quadratic Gaussian technique was also designed and discussed in [21]-[25] for the complex wide area power 
system. 

H-infinity scheme was also discussed in [26]-[27]. These controllers provided good performance 
and robust LFC in the wide area power system. Both SMC, linear quadratic Gaussian (LQG), H-infinity 
schemes worked well under power system with more tuning variables. However, the consideration of 
communication delay has brought a new era for the LFC of MAPS. These delays exist in power systems and 
vary from ten to several hundred on milliseconds. Communication links are used to connect ancillary 
components that are link connecting load frequency sensor point to the remote control unit (RCU) down to 
the generating unit. In reality, the RCU sometimes might experience delay in the receiving ACE signal. In 
practice, time-delay for the controller to receive input ACE signal should not be neglected and this signifies 
the optimal control. Therefore, communication delay should be taking into account in the dynamic modelling 
of complex power system. Newly published paper [28] addressed the LFC for the MAPS with time delay but 
with all system variables are assumed to be measured. In general, not all system parameters can be measured. 
The author work in [29] serves as a motivation to design a real conceptualization controller for the complex 
interconnected power plants with signal delay. In this paper, the LQR combined with Kalman filter technique 
are proposed as LQG method for frequency regulator of the MAPS with communication signal delay. The 
MAPS model is represented in state space with the time-delay to extend the model. This state space model is 
applied to obtain a better gain for the LQG. Since not all values of the LFC model's parameters are 
measurable. Therefore, Kalman filter is utilized for variables estimation. The estimations as well as time- 
delay are input signals. The contribution of the paper is stated below: 

— Communication time delay is considered in the network state space model 

— The proposed controller is designed based on state estimation with feedback gain and time delay 

— Lyapunov theory and linear matrix inequality (LMI) techniques is used to prove the MAPS stability 
— Comparison results are shown between the proposed controller and [29]. 


2. SYSTEM LFC MODEL IN TIME DELAY 

In this section, a conceptualized LFC model taking into details with the time delay is discussed. A 
diagram describing two area power systems is display in Figure 1 [29]. Each area is represented with 
communication delay e*'*. The value x represents the small time delay which might varies from 0.1 to 1 s with 
consideration of two control areas for i,j = 1,2(i # j). Therefore, the modeled system with the information 
of the communication time delays is as (1). 


APrie ij (t) = 20T (Afi(t) — Afj() () 


Other variables deviations are illustrated: 


à 1 Kpi Kpj Kpj 
Afi(t) = — —Afi(t) + F AP mi (t) = T AP rie ij (t) = T AP ai (2) 
pi pi pi pi 
r 1 1 
APqi(t) = — 7, APni(t) + Ta APu(t) (3) 
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AP,i(t) = — ra iO = r, lPi = ioe = 6) r, AO (4) 


where Af;(t) is the frequency deviation that is needed to be controlled, AP,,;(t) is the mechanical power 
change that requires generation-load balance during operation, AP,; is the disturbance affecting the system, 
APrieij(t) is the deviation in the net sum of the power exchange in the communication link connecting both 
area i and j. AE;(t — 1j) which can be denoted as the ACE with communication delays and uj(t) is the 
control input to the system model. 


AE,(t) = KpikKeiAfi(t) + KeiAPtie (t) (5) 
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Figure 1. The structure of the two-area LFC model in a power system with time-delay 


To represent the MAPS, the system state variables Af; (t), AP; (t), AP,;(t), AE;(t), AP rie i; (t) are 
denoted as vector x(t) = [4A,; APtie12 AA,5]" in which AA,, = [Afi (t) AP4í(t) APA(t) AE;(tO] 


and AA,» = [Af2(t) APm2(t) AP,2(t) AE;(t)] The MAPS dynamic equations are derived above and it 
is suggested in the matrix form representation: 


x(t) = Ax(t) + Agi x(t — £1) + Agzx(t — £2) + Bu(t) + 'O(t)y(t) = Cx(t) (6) 


where u(t) = [U1 U2]? is the control signal vector and Q(t) = [4Pg, — APg;]T is the disturbance vector of 
the demands. For l = 1,2 and b = 3,4, the following matrices are defined as: 


Ay A3 04x4 T T T 
: ms p [BE Os]. [f xs 
A=| As 0,4 As |; Aa = [o 0 Tgi 0| ;B = Os B, | r= Ons d; | 
O04 1244 A42 
MS Kpi 0 0 
Tpi Tpi 
0 E K T 
where A, = Ta Ta d ea 0 0 1] ;As- [204 0 0 0] 
1 1 pl 
= o -+ 0 
RiTgi Tgl 
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In the real system, exact measurement of the parameter values is impossible. Therefore, most 
parameters are approximated or estimated. Also, by considering time delay factor, so the dynamic model in 
(6) is further re-written as in (7): 


x(t) = [A + AA(x, t)]x(t) + [Agi + AAgi]x(t — £1) + [Aaz + AA ge |x (t — £2) + 
[B + AB(x, t)Ju(t) + 'Q(t) = Ax(t) + Agi x(t — £1) + Agox(t — £2) + Bu(t) + 
w(x, t)y(t) = Cx(t) + v(t) (7) 


where the unknown matrices AA(x, t) and AB (x, t) denote by system variations, the matrices A and B are the 
nominal values; w(x, t) is called the total variations which can be estimated and we can define it in (8): 


w(x,t) = AA(x, t)x(t) + AAgix(t — &) + AAaax(t — £2) + AB(x,t)u(t) + rA) (8) 


and v(t) is a measurable disturbance. Ideally, some assumptions are stated to describe reality and the 
feasibility of the power system and as well the parameters based load frequency control strategies under 
certain conditions. These assumptions and some lemmas are stated below: 

Assumption 1: If the matrix A, B is controllability, then, AA(x, t), AB(x,t) can be estimated. 

Assumption 2: The lumped w(x,t), which includes delay state matrix and the differential of w(x,t) is 
bounded. i.e., there exist the known scalars y and y such that||w(x, t)|| € y and||w x, t)|| € y, where||. || is 
matrix norm. 

Lemma 1: [13] Let X, Y and F are actual matrices of appropriate dimension with FTF < 1 then, for any 
scalar y > 0, the sequent matrix inequality holds. 


XFY -- YTFTXT < y-1XTX + yy 


3. CONTROLLER DESIGN WITH TIME-DELA Y 
3.1. Linear quadratic regulator (LQR) 

To design the LQG based LFC scheme for the described power system in Figure 1. The derived state 
space model is considered as shown in (7). In practice, the measurement of state variables is difficult and 
expensive because of lack of sensors. Therefore, the Kalman filter technique is realistic to estimate value of 
state variables. In addition, The Kalman filter can reject disturbances in form of white noise v(t) which 
affect the system state. The greater problem is to get a better optimal gain K so that the eigenvalues 
(A — BK) in (9) is places in their desired position and minimized a performance index below. 


J - 5 f[x (00 x(t) + uT (Ru(O]at (9) 


Where Q = E [w(x,t)] having dimension (n x n) and R = E[v(x,t)] having dimension (m x m) both are a 
posited-definite Hermitian or real symmetric matrix. The LQG based controller can be designed in the time 
delay as in (10). 


H(t) = ARCE) + Ag X(t — £1) + Ag; X(t — £2) + 
Bu(t) + Ke) — $(t))$(t) = CEC) (10) 


X(t) is the state variable estimation of the x(t) of the observer and y(t), $(t) are output of the power 
network and observer, respectively. K, is the gain, that is obtained from Riccati equation solution. 


3.2. Proposed controller 
The minimize J in (9) is designed with below controller: 


u = -KÉ(t) (11) 


where K is the optimal feedback gain for the LQR, K = R^1BT9i and 9i is derived also from Riccati 
equation. 


ATR +RA-—RBRIBTR+Q=0 (12) 
Substituting u into (6). 
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x(t) = (A — BK)x(t) + BKX(t) + Aqx(t — £1) + Aa2x (t — €2) + rQ (t) (13) 
Taking the error to be X(t) = x(t) — £(t) and X(t) is calculated as in (14): 
X(t) = x(t) — X(t) = (A — K,C)X(t) + Agi X(t — £1) + AgzX(t — £2) + TA (t) (14) 
where i = 1, 2 are area number and from (1), (2) and (3), we obtained: 


EO] E [4 — BK) 


iol = sepa esc 


(A — K, er Poi Z(t — &) 
^ aie enl íá 


Theorem 1: The system in (1) is asymptotic stable, if there exist matricesR 0, P > 0 and positive scalars u 
and y such that the following LMIS hold: 


(A — BK)'R + R(A - BK) RBK RAs RAg RI 0 0 0 
KTBTR (A—K,C)"P+P(A-K,C) 0 0 0 PAg PAn R 
A aR 0 -u* 0 0 0 0 0 
Tis ò GM ee ee aL) 
0 AT aP 0 0 0 -u 0 0 
0 AT P 0 0 0 0 -u* 0 
0 TTR 0 0 0 0 -u 
Proof: We choose Lyapunov function candidate as shadows 
x(t)]. lee x(t—& il x(t — &) 
veo) - [rol fo sco] * n [ol lo else - mlie- 
t x(t — 8) Qi d £2) 
th lec - E | 0 "1 ES . e d d S 
Then, the time derivative of V (x(t), e(t)) is indicated by (18). 
ses _ [x(t)]' [(A 7 BK)'R + R(A - BK) RBK x(t) 
MAE) = Ba | KTBTR (A — K,C)P + P(A— "e [to 
x(t— &)] [AT aR L] .[x(] [RAg 0 ][x(t— &) 
We = e ie | Beal a E [o 2 Re = &) 
x(t — £) Mook x(t) x(t] [RAg 0 ]1[x(t — &j) 
"ls = 2 ie ar Rl Ral T kol [o 2 [sc — £3) 
x(t) x(t) x(t — &) x(t — &) 
"col E | BA d= [ss 7 A l5 5 [st - &) 
x(t)| [Q: x(t) x(t — €)| [Qi x(t — £2) 
+O) lo. alle 2 - [zc - el lo. o] [st es 
+.:07 (C) T Rx(t) + QT (OT P(t) + x" (RTA) + €T (Prat) (18) 
Applying Lemma 1 to (18), 
X x(t)! | 04 . RBK|[x(t) 
V(x(t), X(t) S Wd lens e 0, | BON + usyi? + yo? + uays? (19) 
hate 
zi. BK)'R + R(A — BK) + uj  RAg AT aR + Mt RAGA’ aR + Ug  RTTTR ; 
(A — K,C)TP + P(A — K,C) + u4 1PAqu AT ai P + u2 1PAgq; AT a2P + u4 1 RTTT Ri; us = ua + uas 
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In addition, the LMIs (16) is equivalent to this inequality: 


0, RBK] _ 


KTB'R 0, =m G0) 


and —Q < 0. Combining (16) and (19), we have: 


V < X25 (Alle) II? nin (21) 
where 0 = [sV17 + ia yz? + uz ys?, the constant value ô and the eigenvalue Amin. Hence, V «0 is succeeded 
with |x| > 2. The system (15) is asymptotically stable. 


4. RESULTS AND DISCUSSION 

In this section, the MAPS with time delay are recommended to test the feasibility of the proposed 
LQG controller. The delays in the area control error signals of £4 and £, are set to be fixed. The parameters of 
the area 1 and 2 generating unit are obtained from paper [29]. 


4.1. Simulation result based on nominal values 

The parameters with their nominal values of the MAPS are used. All state variables are used in 
proposed controller is get from the Kalman filter. At this point, we assume step load disturbances occurring 
on the given system, AP4, = 0.01 pu at 0 s and AP; = 0.015 pu at O s. The results of feedback signals are 
show in Figures 2, and 3 where Figure 2 is frequency deviation of area 1 and area 2, respectively. Figure 3 is 
tie-line power variation between two areas. 

From these figures which show the time variation of the output signal, we can easily observed that 
frequency deviation in both areas converge to zero in 5 s with under shoot are —5 x 107* (pu MW) and 
—2.5 x 1074 (pu MW) in the area 1 and area 2. From these results, it is obvious that proposed controller has 
good performance with soft settling times and small over/under shoots. These results of this controller are 
compared with other method in [29] in next case with the same conditions of disturbance and parameters of 
power system. 


x10? 


Frequency deviation 


5 l l l l l l | l | 
0 1 2 3 4 5 6 7 8 9 10 


Time (sec) 


Figure 2. Frequency deviations (Hz) of the first-area and the second-area 


4.2. Simulation result with random step change load disturbances 

In Figures 4, the signals represent the random step change load turbulences for both area 1 and 
area 2. The random step load turbulences for the two areas network are generated in 200 s for this simulation. 
These load variations are applied to test proposed controller in difference conditions with reference to [29]. 
The simulation results for Af;, 4f; as in [29] are presented in Figure 5 and AP;;, 1? at maximum tolerable 
delay margin in Figure 6, respectively. These results validates that the designed controller gain K is able to 
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achieve stabilization by minimizing the performance index better than [29] at same tolerable delay margin for 


random step load disturbances. 
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Figure 3. Tie line power deviation for power system 
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Figure 4. Load variations 
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CONCLUSION 


As expressed, the LFC based LQR is achieved for the wide area power network with communication 


delays. With regards to lack of sensors, the power system state variables are estimated by the Kalman filter. 
The optimal gain for the proposed controller is selected by solving the Riccati equation so that the control 
parameters are not over tuned. The stability of the system with time delay is again ascertained via Lyapunov 
theory. Results showed better performance to damp transient frequency and ACE than the controller 
compared with it. Also, some of the highlighted limitations from other controllers in literature are solved. 
However, the continuous difficulty in obtaining optimal control feedback gain using some schemes opens a 
new chapter to research area to develop and continuously improving the gain to avoid excessive tuning 
during linearization. 
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